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Abstract 

We consider the propagation of two-photon light in a random medium. We show that the Wigner 
distribution of the two-photon wave function obeys an equation that is analogous to the radiative 
transport equation for classical light. Using this result, we predict that the entanglement of a 
photon pair is destroyed with propagation. 



1 



The propagation of light in disordered media, including clouds, colloidal suspensions and 
biological tissues, is generally considered within the framework of classical optics [1]. How- 
ever, recent experiments have demonstrated the existence of novel effects in multiple light 
scattering, in which the quantized nature of the electromagnetic field is evident. These in- 
clude (i) the transport of quantum noise through random media [2] (ii) the observation of 
spatial correlations in multiply-scattered squeezed light [3, 4] (iii) the measurement of two- 
photon speckle patterns and the observation of non-exponential statistics for two-photon 
correlations [5, 6] and (iv) the finding that interference survives averaging over disorder 
and is manifested as photon correlations, exhibiting both antibunching and anyonic symme- 
try [7, 8]. Thus, there is an interplay between quantum interference and interference due to 
multiple scattering that is of fundamental interest [9-15] and considerable applied impor- 
tance. Indeed, applications to spectroscopy [16], two-photon imaging [17-26] and quantum 
communication [27-29] have been reported. 

In the multiple-scattering regime, the radiative transport equation (RTE) governs the 
propagation of light in random media [1]. The RTE is a conservation law that accounts for 
gains and losses of electromagnetic energy due to scattering and absorption. The physical 
quantity of interest is the specific intensity J(r, k), defined as the intensity at the position r 
in the direction k. The specific intensity obeys the RTE 



which we have written in its stationary form. Here fi a and p, s are the absorption and 
scattering coefficients of the medium and p is the phase function. We note that although 
the RTE is often viewed as phenomenological, it is derivable from the scattering theory of 
electromagnetic waves in a random medium [1, 30, 31]. 

The propagation of two-photon light is generally considered either in free space or, in 
some cases, with account of diffraction [32, 33]. However, understanding the interaction of 
light with matter is central to applications in both imaging and quantum information. In 
this Letter, we consider the propagation of two-photon light in a random medium. We show 
that the averaged Wigner distribution of the two-photon wave function obeys an equation 
that is analogous to the RTE. Using this result, we characterize the loss of entanglement of 
a photon pair upon propagation in a random medium. In this sense, our work builds on the 
well-known duality between partially coherent and partially entangled light, in which loss 
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of entanglement is dual to the gain in coherence with propagation [32, 33]. 

We begin by recalling some important facts about two-photon light. We consider the two- 
photon state and define the second-order coherence function as the normally ordered 
expectation of field operators: 

r( 2 )(r l7 t i; r 2 ,t 2 ) = (^E-^t^E-^^E+^^E+^h)^) , (2) 

where E~ and E + are the negative- and positive-frequency components of the electric-field 
operator with E~ = [E + ]^ . In a material medium with dielectric permittivity £, the field 
operator E + obeys the wave equation [35, 36] 

_ 2 ~, e(r)d 2 E+ 

Here the medium is taken to be nonabsorbing, so that e is purely real. 

The quantity is proportional to the probability of detecting one photon at ri and a 
second photon at r 2 and can be measured in a Hanbury-Brown-Twiss interferometer [37]. 
For the two-photon state \ip), it can be seen that r^ 2 ) factorizes [32, 38] as follows: 

r( 2 )(r l5 tf, r 2 , t 2 ) = M £~(n> h)E-(r 2 , t 2 ) \n) (n\ £+(r 2 , t 2 )E + (v u h) H) (4) 

n 

= (^| E-{r u h)E-(r 2 , t 2 ) |0) (0| E + (r 2 , t 2 )E + (r 1: h) \j>) (5) 
= |$( ri ,t i; r 2 ,t 2 )| 2 . (6) 

Here {\n)} denotes a complete set of states and the two-photon probability amplitude $ is 
defined by 

$(n, t x - r 2 , t 2 ) = (0| E + (r u h)E + {r 2 , t 2 ) \^) . (7) 
Evidently, $ satisfies the pair of wave equations 

which follow from the fact that E + obeys the wave equation (3). We note that (8) is the 
analog the Wolf equations for two-photon light [34] . We will find it convenient to introduce 
the Fourier transform of the probability amplitude $, which is given by 

^(n,^;^,^) = J dM* 2 e^ ltl+ " 2t2) $(ri,*i;r 2 ,t 2 ) . (9) 



Eq. (8) then becomes 

V 2 T fi + k 2 j e(r j )$ = 0, j = 1,2. (10) 

where kj = Uj/c. It is important to note that if $ factorizes into a product of two functions 
which depend upon ri and r 2 separately, then the two-photon state is not entangled. In 
contrast, a fully entangled state is not separable and corresponds to $(ri, oc 5(r\ — r 2 ). 
We now consider the Wigner distribution of $ which is defined by 

W{r, k) = J d 3 r'e lk - r '$(r - r'/2, Wl ; r + r'/2, u 2 ) ■ (11) 

Upon subtracting (10) for j — 1,2 and changing variables according to r x = r — r'/2, r 2 = 
r + r'/2, we see that IV obeys the equation 

\i-V r W+ l -J -0^e-^e(p) [k 2 W(v, k + p/2) - k 2 2 W(r, k - p/2)] = . (12) 

We note that (12) is an exact result which describes the propagation of the Wigner distri- 
bution for two-photon light in a material medium. 

We now proceed to derive the RTE for two-photon light. To this end, we consider a 
statistically homogeneous random medium and assume that the susceptibility rj is a Gaussian 
random field with correlations (r](r)) = 0, (rj(r)r](r')) = C(|r — r'|). Here rj is related to 
the dielectric permittivity by e = 1 + 4irr], C is the two-point correlation function and (• • • ) 
denotes statistical averaging. Let L denote the propagation distance of the field and £ the 
correlation length over which C decays at large distances. We introduce a small parameter 
e = l/(koL) 1 and suppose that the fluctuations in t] are sufficiently weak that C is of 
the order 0(e) and £/L = 0(e). We then rescale the spatial variables according to r x — > 
ri/e, Y2 — > T2/e and define the scaled two-photon probability amplitude <& e ( r i, r 2, ^2) = 
$(ri/e, uji] r 2 /e, oj 2 ), so that (10) becomes 

e 2 V 2 r ^ 6 + fc|$ e = -A7rk 2 V~ev (r,/e) $ e , j = 1, 2 . (13) 

where we have introduced a rescaling of rj to be consistent with the assumption that the 
fluctuations are of size 0(e). If we denote by W e the Wigner distribution of $ e , defined 
according to (11), then (12) becomes 

k • V r W e + ^ (hi - k 2 2 ) W e + = , (14) 

Ze y/e 
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where 

SfW e = 2m J -^e-^/^p) [klW e {r, k + p/2) 

-fc 2 2 W e (r,k-p/2)] . (15) 

We now consider the asymptotics of the Wigner distribution in the homogenization limit e — > 
0. This corresponds to the regime of weak fluctuations. Following standard procedures [30], 
we introduce a two-scale expansion for W e of the form 

W e {r, R, k) = W {r, R, k) + y/eW x (v, R, k) + eW 2 {r, R, k) + • • • , (16) 

where R = r/e is a fast variable. Next, we suppose that /x = (k\ — fc|)/(2fce) = 0(1), which 
corresponds to working in the high-frequency regime. By averaging over the fluctuations on 
the fast scale, it can be seen that (Wo), which we denote by X, obeys the equation 

k • V r X(r, k) + (fi a + /i s )/(r, £) = n a f d 2 k'f(k, k')X(r, k') . (17) 

Here the absorption coefficient fi a , scattering coefficient /x s and scattering kernel / are defined 
by 

fji a = (kf - jfc|) J C{k{k - k'))d 2 k' + in , (18) 
I!, = k\(k\ + k 2 2 ) J C{k{k - k!))d 2 k' , (19) 

where f(k, k!) = C(k(k— k')) is normalized so that J /(k, k')d 2 k' = 1 for all k. We note that 
this normalization is consistent with the statistical homogeneity of the random medium since 
C(k(k — k )) depends only upon the quantity k ■ k'. We will refer to (17) as the two-photon 
RTE. Evidently, (17) is the analog of the classical RTE. However, the physical interpretation 
of (17) requires some care. In contrast to the specific intensity, the quantity X is not real- 
valued and is not directly measurable. Nevertheless, by inversion of the Fourier transform 
(11), we find that X is related to the average two-photon probability amplitude by means of 
the formula 

(5(n,r 2 )) = /(gje-<-.->x(^,k) , (20) 

where the dependence on the frequencies u\ and u 2 has not been indicated. 

We now explore some physical consequences of the two-photon RTE. In particular, we 
examine the propagation of entanglement. We begin with the case of a deterministic medium 



in which the permittivity e is constant. We consider the half-space z > and assume that 
the two-photon Wigner distribution Iq is specified on the disk of radius a in the plane z = 
in the ingoing direction. That is, 

A5(k — kn) if k • z > and \p\ < a 
X (r,k)=<| (21) 

otherwise , 

where A is constant and p is the transverse coordinate in the z = plane. Making use 
of (20), it is readily seen that 

( 2nk$A si f° lpi - p f ] if \ Pl2 \<a 
$( Pl ,0;p 2 ,0)= ° ^p*\ l^ 2 '- (22) 

I otherwise , 

which corresponds to a transversely entangled two-photon state. To propagate X into the 
z > half-space, we make use of the formula 

X(r, k) = f d 2 k' [ d 2 r' z ■ k'G(r, k; r', k')X (r', k') . (23) 

J Jz'=0 

Here G is the Green's function for the two-photon RTE (12), which is given by 

G(r, k, r', k') = y^A* ~ &)* ( fc " j^f ) ■ (24) 

We can now compute the two-photon probability amplitude $. For simplicity, we assume 
that k\ — ki — ko and that the points of observation ri and r 2 are on-axis, with ri = r 2 = 
(0, z). Carrying out the integrations in (23) and making use of (20) we find that 

5(o, *M = 2 *»- (s) ■ < 25 > 

It can be seen that the entanglement of the photon pair is destroyed with propagation. We 
note that the diagonal part of the coherence function r^ 2 ^(r, r) = |$(r, r)| 2 is proportional 
to the probability of two-photon absorption at the point r. Next, we consider the case of 
a random medium. At first, we will make use of the diffusion approximation (DA) to the 
RTE, which is widely used in applications. The DA neglects the angular dependence of 
the Green's function. It holds in the limit of strong scattering and at large distances from 
the source. Within the accuracy of the DA, the Green's function for the RTE is given by 



G(r, r') = exp(-fi:|r - r'\)/{^D\r - r'|). Here « = ^JD, D = l/[3(/i a + (1 - g)fi s )} and 




FIG. 1. (Color Online) Dependence of (<5) on the distance of propagation z for different off-axis 
distances p. The scattering kernel / was taken to be of the Henyey-Greenstein form with g = 0.9, 
Ha = and p, s = lOOfco- 

g — J It ■ k'/(k, k.')d 2 k'. Carrying out the integrations in (20) and (23), we find tha the 
average two-photon probability amplitude is given by 

< 8 ( r ^>>= ^V iD f°' r '~i r2 ' ) < 26 > 

2U(2tt) z \Ti — r 2 | 



x 



[°° ^_ J x {qa)JMPx + p 2 |/2)e-v^(-+-)/ 2 , 
Jo \/q z + K 



\/q 2 + k 2 

where r = (p, z). In the on-axis configuration, we find that 

{m Z] °' " )} = 2D^Y 2 [^+^ ~ z \ ■ ( 27 ) 
As above, the entanglement of the photon pair is destroyed with propagation. Finally, we 
consider the propagation of entanglement in the transport regime. The Green's function 
for the two-photon RTE can be obtained using the method of rotated reference frames [39- 
41]. The calculation of the average two-photon probability amplitude is presented in the 
supplementary material. In Fig. 1 we plot the ^-dependence of ($(r, r)) for various values 
of the off-axis distance p with k\ = &2 = ho- Once again, the entanglement of the photon 
pair is lost with propagation. 

We close with a few remarks, (i) It is possible to derive the analog of the RTE for single 
photons. Not surprisingly, this equation has the form of the classical RTE (1). We plan to 



present (ii) Although in our model the electromagnetic field is quantized, the interaction of 
the field with the scattering medium is treated classically. It would be of interest to extend 
our results to the case in which the medium consists of a collection of two- or three-level 
atoms. In this manner, it should (in principle) be possible to understand the transfer of 
entanglement from the field to the medium [42]. Evidently, the calculations that we have 
presented do not account for this effect, since we have taken a macroscopic approach to the 
quantization of the field [35, 36]. (iii) Finally, applications to imaging and communication 
theory may be envisioned. In the former case, there has been extensive use of the classical 
RTE for imaging in random media. It may be anticipated that experiments with two-photon 
light may enjoy some advantages, as has been suggested for the case of quantum optical 
coherence tomography [43, 44]. In the latter case, there has been considerable interest in 
the use of quantum states of light for communication [27, 29, 45]. It would be of interest 
to understand the effect of a random medium, such as the atmosphere, on the capacity of 
quantum information systems [28]. 
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